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SIMULATOR FOR STATIC LIQUID CONFIGURATION IN PROPELLANT 
TANKS SUBJECT TO REDUCED GRAVITY 
by William A. Olsen 
Lewis Research Center 


SUMMARY 

A simple device has Been investigated and constructed which can approxi- 
mately simulate (show) the isothermal equilibrium liquid configuration in pro- 
pellant tanks of various shapes and sizes that are subject to reduced gravity 
ranging from 1 to near "zero” gravity. Experimental results compare favorably 
with drop- tower results. 

The simulator, an application of the "Hele-Shaw Cell, " is basically com- 
posed of two parallel, closely spaced, flat glass plates that are separated by 
rubber gaskets, which are cut to form the two-dimensional tank boundary. A 
given amount of liquid is placed within these boundaries. The angle of the 
plates to the horizontal is varied to obtain various effective gravity fields. 
The basic principle of the simulator can be roughly stated in the following 
manner: When two parallel, flat plates are closely spaced, the two-dimensional 

liquid shape in the plane of the plates is only affected by the gravity compo- 
nent and the tank boundaries in that plane. Thus, when the plates are perfectly 
horizontal, the liquid should assume the zero-gravity shape 1 for the particular 
two-dimensional tank boundaries and liquid used between the plates. 

The design and operating criteria for the device that allow it to simulate 
the static liquid shape in two-dimensional tanks is obtained analytically and 
substantiated by experiments. The ability of the device to simulate three- 
dimensional tanks approximately is indicated experimentally. Sample photographs 
were taken of the equilibrium liquid configuration within the simulator for var- 
ious gravity fields, tank sizes, and configurations. 


INTRODUCTION 

The NASA Lewis Research Center is conducting studies of the behavior of 
propellants in storage tanks while exposed to reduced gravity (ref. l). This 
condition will be encountered by space vehicles during planetary orbital mis- 
sions and during periods of coast on interplanetary space missions (ref. 2). A 
knowledge of the liquid configuration will be needed to solve problems related 
to effective tank venting, liquid transfer, orientation control of the vehicles, 

^or a given tank and liquid, zero gravity is that gravity below which 
changes in gravity have a negligible effect on liquid configuration (i.e.. Bond 
number < 10“^ ) . 


and heat transfer (ref. l). 

Reduced gravity studies of the static liquid configuration in propellant 
tanks have been generally limited to specific analytical studies (refs. 3 
and 4) and experiments (refs. 5 and 6). A device that quickly, simply, and 
inexpensively shows the approximate shape of a liquid in a tank subject to re- 
duced gravity could greatly help the engineer and designer in this new field. 
Many different experimental methods have been used to obtain the reduced- gravity 
liquid configuration, such as drop towers, zero-gravity trajectories in air- 
craft, and a simulator using the same density but different surface tension liq- 
uids. In this investigation another approach was taken that involved a two- 
dimensional simulator to determine the approximate static liquid configuration 
within propellant tanks of varying shapes and sizes, which were subject to ef- 
fective gravity fields from near zero gravity to 1 g. 

This report describes this simulator and shows some of the results obtained. 
The simulator is composed of two parallel, closely spaced, flat glass plates, 
separated by rubber gaskets, which form the two-dimensional tank boundary. Liq- 
uid is placed within these boundaries. The angle of the plates to the horizon- 
tal is varied to obtain various effective gravity fields. 

The idea of using two parallel glass plates for fluid studies is not new. 
Hele-Shaw (ref. 7) devised a similar arrangement to study two-dimensional po- 
tential flow around objects. Reference 8 later used this same method to simu- 
late the stability of the interface between two different liquids as they flowed 
through a porous media. In this reference, it is suggested that the effect of 
gravity would disappear when the plates were horizontal. 

The basis for the device in this application is derived and discussed. 

Some liquid configuration data are obtained for static situations. 


ANALYTICAL BASIS FOR SIMULATION OF STATIC LIQUID SHAPE IN TANKS 

In this section the design and operating conditions for a two-dimensional 
simulator of the static liquid shape in containers of various sizes subject to 
varying gravity fields will be derived. Theory for scaling simulator results to 
two- and three-dimensional containers of various sizes in various gravity fields 
will be discussed. The analysis that follows is basically an application of 
the type of analysis found in a number of references, for example references 3 
and 4. 


General Equations 

Consider the energy of a fixed amount of liquid at rest, at constant tem- 
perature, with constant properties, In a simply formed tank where surface ten- 
sion forces are significant. The shape of the liquid depends on the sum of the 
surface free energies and the potential energy of the liquid (ref. 9) 
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■where 

p e = P Z - p f 

(All symbols are defined in appendix A. ) 

Expand equation (l) to include the surface free energies, respectively, of 
the liquid- solid (LS), solid- vapor (SV), and liquid- vapor (LV) interfaces of 
figure 1 to obtain 

1 ' °LS ff dS + °LV ff as + » sv ff dS 4- g 0 p e f ff nz dV 
b LS S LV b SV v 

(2) 

The Dupre-Young equation (refs. 3 and 10 ), which relates the surface free 
energies per unit area to the contact angle 2 0, is given as follows: 

°SV ~ a LS = a LV COS 6 

This equation strictly holds for 0° < 0 < 180°, hut it is assumed to he suffi- 
ciently accurate for this analysis even when 0=0. 

Comhine equation (2) with (3) to eliminate a^. Then use equation (4), 
which breaks up the area of the container: 

Sb ‘ff as + ff ds <*> 

S SV S LS 

These substitutions result in 

1 - a L¥ // dS + °LV COS 9 ff dS + °LS^b + So p e f ff nz dV 
S LV S SV v 

(5) 

The configuration of a fixed amount of liquid in a tank is in stable equilibrium 
when this energy I is a minimum (ref. 12). Mathematically, this requires the 
variation of I to vanish (i.e., 81 = 0) subject to the constraint of constant 
liquid volume. Performing the required variation operation upon equation (5) 


^Available data from isothermal drop-tower experiments with no flow appear 
to indicate that the contact angle 9 is a constant (ref. 11 ). 
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results in the following equation (eq. (6a) ), which together with the con- 
straint equation (eq. (7)) describes the liquid configuration in a given three- 
dimensional tank subject to a given gravity field 





(6a) 


and 


dV = constant = V (7) 

It is possible at this point to simplify equation (6a) by eliminating the term 
a Ls£>b in an actual three-dimensional case, .since it is a constant whose vari- 
ation vanishes. This equation would then become 
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(6b) 


However, this term must be kept for the subsequent analysis of the simulator be- 
cause it may not be constant in that case. 


Two-Dimensional Simulator 


Equations (6a) and (7) will now be applied to the case of a two- 
dimensional tank model in a 1-g field which is bounded top and bottom by glass 
plates and around the edges by a two-dimensional model of a tank configuration 
(see fig. 2). In this situation, the surface energy is made up of the surface 
energies due to the two-dimensional model boundaries, the gas, and the top and 
bottom plates. Assume the surface tension properties are the same for edges and 
plates, then the equations governing the liquid configuration in this case are: 


J LV 


ffs 

S LV 


dS + a - ky cos 6 


iff AS + 2 f f as ) 

\ J •'Ssvb J • / 3 S¥p / 

+ + 2S p) + SoPe ffl nz dV 


= 0 ( 8 ) 


where 


8b = S B + 2S p 


and 
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Equation (8) is subject to the constraint 


"SVP 
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dV = constant = V 


( 9 ) 


By analogy with the equations for a three-dimensional tank (eqs. (6a) and (7)) 
the equations for a two-dimensional tank would he 


A £ 

' U LV 


ds + Ojjj cos 8 


ds + SoPe m dA^ 


= 0 


( 10 ) 


and 




dA = constant = A 


(id 


The desired result is to have the simulator, a three-dimensional device, simu- 
late the liquid shape in two-dimensional containers^. The question then is what 
assumptions are necessary to transform equations (8) and (9) into the form of 
euqations (10) and (ll). These assumptions are given as follows: 

(1) The plate spacing is small compared with the local characteristic width 
of the interface in the model. Usually this characteristic width is taken as 
the tank or capillary tube diameter. 

(2) The plate spacing is less than the maximum spacing required for the 
liquid to contact both plates (see appendix B). 

(3) The plates are flat and parallel (i.e., £ = const.). 


By these assumptions the liquid area A(x, a) (fig. 2), in the planes paral- 
lel to the plates is essentially unchanging in the direction normal to the 


plates (i. 
of volume 


e *> f ds and ff dA 
dV can be taken as 


are not functions of (3). Also, the element 


dV = | dx da = £ dA 

Consider now the potential energy term of equation (8). Consider the 


A two-dimensional container is formed by infinite translation of the two- 
dimensional (plane) shape in the direction normal to the plane of the two- 
dimensional shape. 
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parallel plates to "be inclined at an angle cp to the horizontal (see fig. 2). 
Then by using the preceding assumptions this term becomes 


^ SSL nz dV = g Q p e ffJL na sin 9 dV = g Q p e | a ma dA 

where m = n sin cp, z = a sin cp, and dV = £ cbc da. 

By the previous assumptions the integrals can each be integrated once in 
the p direction. Integrating each term of equations (8) and (9) results in 


V = constant 




H 


dA - constant = A 


and similarly 


// 

sv 


dA = Ag^. = constant 


°LV If dS = a LV f ds y dp = a LV | j) ds 

S LV C LV 0 C LV 

2( a LV cos 0) If dS = 2(agy cos 0) // dA = 2( cos £)Agy = constant 

A sv 


"SVP 


Oj-y COS 9 


a 


dS = cos 0 


SVB 


f as / 

°SV 0 


dp = cTjpy cos 9 1 (f ds 


°sv 


Finally, 


■W s b + 2 V - “ls'V + ZA r> 


Substitute these relations into equations (8) and (9) and divide by | to ob- 
tain 


LS 


(% + ~r) + 20 LV C0S * -T * °LV f as 

U LV 

/ ds + g n P P ® ® ma dA 

•/ 


CTjy cos 9 


=0 ( 12 ) 
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The first two terms in equation (12), not involving integrals, are con- 
stants for a given case and the variation of a constant is zero; therefore, the 
equation can he reduced to 
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ds + 


°LV 


cos 0 




Equation (13), together with the constraint 


(13) 



dA = constant 


= A 


(14) 


describes the shape of the liquid in the simulator. The equations describing 
the liquid shape in the simulator (eqs. (13) and (14)) and for a two-dimensional 
container (eqs. (10) and (11)) are of the same form if the assumptions used in 
the derivations are correct. By inspection, the two-dimensional liquid shape in 
the simulator is only affected by tank boundaries and surface energies in the 
plane of plates (a, x plane) with the magnitude of gravity taken as a component 
in that plane (m = sin cp). The two-dimensional liquid shape is not affected by 
the plates. This result assumed the simulator to be designed with flat, paral- 
lel plates that are separated by a distance that is less than the maximum re- 
quired for wetting and also required the characteristic interface dimension to 
be much greater than the plate spacing. 


In order to make the liquid configuration in the simulator and in a larger 
two-dimensional container the same, the terms in the governing equations must be 
of the same magnitude, term by term (i.e., scaled), in addition to being of the 
same form. Scaling equations for magnitude is best handled by nondimensional- 
izing the variational equations describing the two-dimensional liquid configura- 
tion (eqs. (13) and (14)). 


Equations describing the two-dimensional simulator and a two-dimensional 
container (eqs. (13) and (14)) are nondimensionalized by the following substi- 
tions 


a = ?£ s 
a = a,zl 


S LV " C LV £ s 
S SV = ^SV^s 


£ 

LV 


do + cos(0 ) 


and 


£ 


dO + 


SV 



da 


= o 


(15) 



(16) 
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where 


m = n for two-dimensional case 

m = sin cp for simulator and n = 1 

The nondimensionalization has "been done so that the nondimens ional variable in- 
tegrals (lower-case script letters) vary within a unit order of magnitude, while 
the constant coefficients (nondimensional groups) in front of the integral 
signs, containing characteristic quantities, represent the magnitude of the 
term of which they are a part. It was assumed that adequate scaling could be 
achieved if the characteristic quantities were taken as a characteristic dimen- 
sion of the tank, such as the tank diameter. The validity of this nondimension- 
alization must be checked by an experiment where the container size, configura- 
tion, liquid properties, and gravity field are varied. The term in brackets in 
equation (15) is commonly called the Bond number B. It represents the ratio 
of gravitational to surface energy. Examination of equation (15) reveals that 
the gravitational energy term (third term) will not appreciably affect the liq- 
uid configuration compared with surface energy terms for small Bond numbers 
(B « 1), while the surface energy terms (first two terms) have little affect 
for large Bond numbers (B » l). The former case is commonly referred to as 
"zero" gravity. The latter case is common to our everyday experience (e.g., a 
glass of water on a desk). It should be expected that the important effects on 
the liquid configuration of varying the Bond number can be determined over a 
Bond number range of 10^ > B > 10“^. 

According to equations (15) and (16) the two-dimensional simulator should 
be able to show the shape of the liquid in a different sized two-dimensional 
container if the following scaling conditions are met: 

(1) The contact angles 0 are the same for each two-dimensional container. 

(2) The Bond numbers are of equal magnitude for each. 

(3) The containers and liquid content are geometrically similar (ref. 13); 
that is, the corresponding two-dimensional tank dimensions are proportional and 
the liquid content fraction, based on area, is the same. 

Numerous authors (e.g., refs. 1 and 3) have noted these same scaling re- 
quirements for three-dimensional containers. 


Simulation of Three-Dimensional Tanks 


In this section an attempt will be made to determine the three-dimensional 
tanks that can be studied in the two-dimensional simulator. Nondimensionalize 
equations (6a) and (7) by the following substitution: 


= f £ a 

II 

CO 

^LV^a 

il 

CM 

s sv - 

^SV £ a 
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The equations describing the actual case become 


and 




( 17 ) 



(18) 


It has been common practice for engineers to solve two- dimensional axisym- 
metric problems by assuming them to be approximated by a two-dimensional plane 
problem so that more simple analytical or graphical solution methods could be 
employed. This assumption is also made here, with the knowledge that the ap- 
proximation is best away from the axis of rotation. In other words the assump- 
tion is made that the liquid configuration in the two-dimensional simulator is 
approximately the same as the two-dimensional projection (fig. l) of the three- 
dimensional liquid configuration in a geometrically similar axi symmetric three- 
dimensional tank. Based on the previous scaling results, it would be expected 
that the actual three-dimensional case and the simulator case should be scaled 
if 


(1) The contact angles are equal 

(2) The Bond numbers are equal 


= 0 s 


B s = B a 


or 




(19) 


( 20 ) 


In addition, the two-dimensional projection of the three-dimensional tank 
should he geometrically similar to the simulator two-dimensional tank and the 
liquid content fractions, based on two-dimensional projected area or three- 
dimensional volume (simulator shape rotated 360°), should also he geometrically 
similar. In order to simulate a three-dimensional tank containing a constant 
percent full hy volume, the amount of liquid in the two-dimensional simulator 
must he varied some as the effective gravity field m varies. 

The validity of the assumptions made in this section must he checked hy 
experiment. Such an experiment would compare drop-tower static liquid configu- 
ration data for three-dimensional models (surfaces of revolution) to correspond- 
ing simulator data where the effective gravity field, liquid fraction, and 
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container size and configuration are varied. 


Summary of Analysis 

The analysis indicated that the two-dimensional shape of the liquid in the 
plane of the plates of the properly designed two-dimensional simulator should 
only be affected by the magnitude of gravity, taken as a component, and tank 
boundaries in that plane. A properly designed simulator requires two parallel 
flat glass plates that are closely spaced so that the liquid contacts them both 
and the spacing is far smaller than the local characteristic width of the inter- 
face. The analysis indicated that the properly designed simulator should simu- 
late (show) the shape of the liquid in a two-dimensional container. In addition 
it was assumed that the projected two-dimensional liquid shape in three- 
dimensional containers (surfaces of revolution) may be approximately simulated 
by the two-dimensional simulator. Experiments must be performed to validate the 
assumptions of the analysis. 


APPARATUS AND PROCEDURE 
Apparatus 

Based on the theoretical results previously described, a small simulator 
was constructed (fig. 3) whose final design and operational procedure, described 
subsequently, resulted from a series of tests for this purpose. Basically, the 
device consists of two parallel, flat, glass plates separated by shaped rubber 
gaskets, which also act as a seal. The gasket is cut to shape to form a scale 
model of the two-dimensional boundaries of the tank shape of interest. 

The glass plate and rubber "sandwich" is clamped to the flat-bed plate 
(fig. 3) that can be pivoted about an axis. Ideally, when the bed is perfectly 
level the liquid in the test section (the space between the glass plates and 
within the gasket boundaries) should assume the zero-gravity shape. “When the 
bed is at an angle to the horizontal cp accomplished by placing a machined 
spacer under the bed leg, this liquid should assume the corresponding shape for 
a fraction of 1 g. The bottom glass plate is clamped to the bed by a vacuum be- 
tween them. The upper glass plate is clamped by another glass-plate gasket- 
sandwich arrangement, which tends to decrease clamping errors. Closely machined 
spacers (accuracy ±2.5x10*"^ cm) control the spacing | between the glass plates 
and keep the plates nearly parallel. Liquid contact with the plates is assured 
by choosing a plate spacing £ that is less than the "break- away" spacing, 
which was determined experimentally to be l max = °*28 centimeter for ethyl al- 
cohol (see appendix B). All the data in this report were taken with 
£ = 0. 152- centimeter spacers. The glass plates are flat to 2.54x10-5 centimeter. 
The device is leveled by three fine-thread leveling screws. 

Ethyl alcohol (200-proof ethanol) has been used as the test liquid for most 
of the experiments because it is easy to handle and photograph. The alcohol is 
placed in the test section by a hypodermic needle that pierces the rubber gasket. 
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Another hypodermic needle is used as a vent. For this experiment the simulator 
was housed in a room that could not he considered to he a clean room; it was 
qualitatively about the level of cleanliness of an average household. The 
plates were cleaned with alcohol and a ruhher squeegee every time there was a 
change in the model. 

Ethyl alcohol has a contact angle of 0=0° on glass and ruhher. Nonwet- 
ting liquids could, in principle, he used in the simulator. Mercury (0 ~ 125°) 
was tried; however, dust particles and/or oxidation appeared to affect the flow 
of this liquid detrimentally to such a degree that this phase of the experiment 
was dropped. The alcohol exhibited no such difficulties. 


Procedure 

The two-dimensional tank model centerline is oriented with respect to the 
gravity vector (downhill direction) to he consistent with the actual case to he 
simulated (fig. 2). The simulator is leveled (the leveling procedure is dis- 
cussed in the next section) and an amount of liquid, which will give the geomet- 
rically similar liquid content as found in the actual case to he simulated, is 
injected into the test section. The bed is then stopped from the "leveled" 
position to a desired bed angle cp, which is between 10"^ to rt/2 radians. 
Stopped positions are obtained by placing spacers between the base plate and the 
bed leg. The equilibrium (static) liquid shape for various bed angles was ob- 
served and photographed. 


Leveling Procedure 

It is desirable for the test section to be level to a high degree of accu- 
racy (cp = 10“ ^ to 10“^ rad), and the plates parallel to a comparable accuracy, 
for low-gravity-field simulation; therefore, a special leveling procedure was 
developed. Leveling is accomplished by adjusting the leveling screws (fig. 3) 
until an air bubble injected into an unobstructed liquid area of the test 
section, moves unmolested by dust, etc. in the opposite direction from the de- 
sired gravity vector and at a very slow speed. Prom the bubble speed U the 
bed angle can be determined from figure 4 which was derived in appendix C. For 
example, a bed angle of cp = 10“ 4 radian will result in a bubble velocity of 
U ^ 0.3 centimeter per minute in ethyl alcohol for a spacing of | = 0.152 cen- 
timeter and a bubble radius of a = 0.64 centimeter. To aid in this leveling 
procedure, the bottom glass plate can be initially leveled with a high accu- 
racy level (cp = 0.835xl0“ 5 rad). It has been experimentally determined that 
bubble leveling is required, in addition to leveling by a level, at a small bed 
angle because this method automatically compensates for nonparallel plates 
caused by irregularities and also because the effective bed angle can be ade- 
quately estimated by observing the bubble speed and direction. In other words, 
the test section itself is used as a level. Machining inaccuracies of the 
clamps and spacers and the presence of dust particles limited the minimum meas- 
urable bed angle to about 5xl0“5 radian for ethyl alcohol. Once the bed is 
leveled by this method the experimenter, can then set other bed angles (i.e., 
stopped positions) from about 10"^ to jt/2 radians with good accuracy by merely 
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inserting appropriate machined spacers between the base plate and the bed leg. 
After a reduced gravity position or angle is set, the bed can be pivoted to 
another angle (or gravity field) and returned to the reduced gravity position 
repeatedly. 


Operational Capabilities 

The operational capabilities of the two-dimensional simulator discussed 
herein are summarized as follows: 


Time duration possible for experiments Limitless 

Measurable bed angle range, rad re/ 2 > cp > 5xl0~5 

Range of model size (diam), cm 25 > D s > 2.5 

Contact angle (with air), deg 0 (ethyl alcohol) 

Bond number range 10 4 > B s > 2X10“^ 


RESULTS AND DISCUSSION 
Scaling 

The theoretical analysis indicated that the equilibrium liquid shape in a 
properly designed two-dimensional simulator should agree with the liquid shape 
in another geometrically similar two-dimensional container of different size if 
the contact angles and Bond numbers for each case were equal and the liquid con- 
tent was geometrically similar. Experiments were performed to determine the 
validity of this theoretical result. 

To check the Bond number scaling over a range of gravity fields, the simu- 
lator itself was used. Two different sized, geometrically similar, circular 
tanks, half full of liquid, were subjected to varying bed angles (effective 
gravity fields) in the simulator. Figures 5(a) to (c) compare the shape of the 
liquid in a 7.6-centimeter-diameter circle to the shape occurring in a 12.7- 
centimeter- diameter circle at the same bed angle (effective gravity field) and 
also at a bed angle required for correct Bond number scaling. The dark circle 
seen in the photographs is a reflection of the out-of-focus camera lens, and the 
two horizontal lines are channels in the bed plate for vacuum distribution. The 
liquid region of the photographs has been darkened to help the reader distin- 
quish the liquid from the gas. The darkest curved line is the interface while 
the neighboring lighter curved line is a reflection of the darker line. The ef- 
fective gravity field is directed toward the bottom of the page. Figures 5(d) 
to (f) show the result of a similar comparison when the circular tank contains 
a channel. Figure 5 demonstrates that the shape of the liquid in the 7.6- and 
12.7-centimeter-diameter circles is the same only when the Bond numbers are 
equal. This result is shown most clearly in the center columns of this figure, 
while the remaining photographs must be superimposed to see this clearly. The 
practical difficulty with this method of comparison is a result of the fact that 
changes in the Bond number will greatly affect the liquid shape only at the in- 
termediate Bond numbers, such as those of the center column, where the surface 
forces and gravity forces are of the same order of magnitude. 
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Based on figure 5, it can be stated that the Bond number scaling is appar- 
ently a valid method to scale the simulator results to a different sized geo- 
metrically similar two-dimensional tanks subjected to varying gravity fields. 
Because only one liquid (ethyl alcohol) has been used in the simulator, contact 
angle and surface tension were not varied so that a complete check of the scal- 
ing parameters for the simulator, indicated in the analysis, could not be per- 
formed. Bond number contact-angle scaling is generally accepted (refs. 1 and 3) 
for three-dimensional containers so that this incompleteness is probably not im- 
portant. 


Simulation of Three-Dimensional Tanks 

In this section photographs of several of the two-dimensional static liquid 
configurations that occurred in the simulator are compared with the two- 
dimensional photographs of the liquid configuration that occurred in correspond- 
ing geometrically similar three-dimensional surface of revolution tanks in the 
drop tower (ref. 11). Available drop-tower data are presently limited to very 
small Bond numbers (B a < 10“ 3 , where the gravity field is about 10“ 6 g and the 
containers are small). Liquid configuration data at the other extreme of very 
large Bond numbers (B a > 10^, e.g., where the gravity field is 1 g) requires 
only desk top observations. As a consequence, data from the two-dimensional 
simulator and three-dimensional drop-tower liquid configuration (two-dimensional 
photographs) can only be compared at these two extremes and only for those tank- 
shapes that have been investigated In drop tower studies. 

Figures 6 to 8 compare the liquid configuration in three different tank 
configurations at various liquid content fractions (based on volume) and at zero 
gravity ( very low Bond number ) . 

The two-dimensional containers in the simulator are filled to a level so 
that the rotation of the two-dimensional area about the centerline will result 
in the same fractional content (percent full by volume) as for the three- 
dimensional drop- tower case at a given Bond number. This would normally be done 
iteratively by first considering the percent full by area equal to the percent 
full by volume and then tracing the resulting two-dimensional liquid shape that 
occurs at that Bond number. From this sketch, the rotated percent full by vol- 
ume could be calculated. By continually correcting the amount of liquid in the 
simulator until the desired percent full by volume is attained, by rotation of 
the two-dimensional area, it is possible to obtain a geometrically similar con- 
tent of liquid (by volume) in the simulator. 

Figure 6 indicates good agreement between the simulator and the drop tower 
considering the fact that the air bubble can be anywhere at zero gravity and the 
Bond number for the simulator data was not quite as low as for the drop-tower 
tests. There is good agreement between simulator and drop-tower data for the 
cylinder and cone (figs. 7 and 8) except for some liquid that formed at the top 
of the cylinder and cone for the simulator equilibrium data. This discrepancy 
in equilibrium shapes is probably a result of some "edge effect" for the simu- 
lator. In the analysis it was essentially assumed that the plate spacing was 
much smaller than the local interface width. This assumption would fail where 
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a narrow width of liquid occurs in the simulator such as the thin layer that 
slowly flows up the walls of the cylinder and cone in the simulator and forms 
liquid at the top of these models. The flow up the walls is very slow so that 
the time to reach the simulator equilibrium shape is considerably longer than 
the time required to reach the pseudo-equilibrium shape (before flow up the 
walls) of figures 7 and 8. The size of the pseudo-equilibrium shape then grows 
smaller as a result of this flow until the simulator equilibrium shape is at- 
tained. In most cases this edge effect should prove to be no problem because 
the flow is very slow in the narrow layers so that the pseudo-equilibrium shape 
should be recognized readily. No formal comparison is required at very high 
Bond numbers (B a > 10 4 ) since the liquid configuration in three-dimensional con- 
tainers is well known to the reader (the interface is essentially flat, perpen- 
dicular to the gravity vector). This flat shape also occurs at correspondingly 
high Bond numbers in the simulator. 


Equilibrium Liquid Configuration Data 

Photographs of the equilibrium ( static") shape of totally wetting liquid 
(ethyl alcohol, 0 = 0°) in the simulator, for various two-dimensional tank con- 
figurations and bed angles, were taken. These photographs (figs. 9 to 15), 
when scaled, can be taken to portray (show) approximately the two-dimensional 
projected liquid shape (fig. l) in corresponding three-dimensional tanks (sur- 
faces of revolution) when the edge effect contribution is ignored. Each figure 
is for a constant projected liquid area fraction rather than for a constant liq- 
uid volume fraction. 

Figures 9 to 11 show the approximate two-dimensional simulator equilibrium 
shape of the liquid (shaded areas) at various Bond numbers (based on diameter) 
in a 12.7-centimeter-diameter circle (sphere) and a rectangle (cylinder, 

L/D = 1.55 and 9.5 cm in diam). The lowest Bond number results were necessarily 
obtained with a 2. 54- centimeter model because the bed angle could not be made 
small enough to reach the desired Bond number without reducing the model size. 
These two-dimensional results can be scaled by the Bond number scaling parameter 
to larger geometrically similar two-dimensional (three dimensional) tanks at 
various gravity fields. For example, figures 9 and 10 show the liquid shape 
approximately in a 2. 5-meter- diameter circle (sphere), partly filled with ethyl 
alcohol, at gravity fields of 2x10" 4, 2. 7x10" 6, 2.3xl0"6, 2. 7x10" 7, and 1x10 "8 
in the order of decreasing Bond number. Figure 12 depicts the effect on the 
liquid configuration of the diameter N of an open-end channel (cylinder) inside 
a rectangular (cylindrical) tank. The liquid tends to fill the narrower region. 
Figures 11 and 12, which are simulator equilibrium data, show further examples 
of the result of the edge effect problem previously discussed. Actually the 
edge effect problem occurs in all the figures (i.e., the boundaries are wetted) 
but it does not show up in the photographs except at corners where liquid can 
build up. Figure 13 demonstrates the equilibrium shape of the liquid in a cir- 
cle (sphere), with an open-end channel (cylinder), subject to varying gravity 
fields. Note that the liquid shape tends to become independent of the gravity 
field direction as the Bond number becomes very small. The results of figures 
9 to 12 are in agreement with those of reference 1 for similar situations. 
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The ability of the two-dimensional simulator to simulate the liquid shape 
in three-dimensional tanks is demonstrated in figure 13. The results of the 
smallest hed angle are compared with the drop tower results in figure 14. Fig- 
ure 15 depicts the shape of liquid at various Bond numbers in a typical rocket 
propellant tank that is 14 percent full (based on projected area). 


CONCLUDING REMARKS 

The two-dimensional simulator discussed in this report appears to be able 
to show the approximate equilibrium liquid shape in many surf ace-of-revolution 
containers over a range of tank size and gravity field when scaled by the Bond 
number and the proper procedure is used in its operation and interpretation of 
its data. This device has some advantages over other methods when the approxi- 
mate liquid shape over the significant Bond number range is desired. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, October 7, 1965. 
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APPENDIX A 


SYMBOLS 

A liquid contact area with either upper or lower plates, sq cm 
A,p area of plate within model Boundaries, sq cm 
a Bubble radius, cm 

d nondimensional area, = A/£^ 

B Bond number 

C length of curve, cm 

effective drag coefficient for bubble 

0 nondimensional length variable, = s/Z 

D characteristic tank diameter, cm 

e constant 

F-p. bouyant force, dynes 

-D 

drag force, dynes 

g Q constant acceleration due to gravity, 981 cm/sec 2 

1 energy, dyne -cm 

K constant 

£ characteristic length, cm 

M, N baffle spacing defined in fig. 12 

m n sin cp, number of gravities component in opposite direction of a 
(fig. 2); n = 1 for simulator, therefore, m = sin cp 

n number of gravities in opposite direction of z (fig. 2) 

modified Reynolds number, = (U|/v^ )( £/a)^/ 2 

r radial coordinate 

S surface area, sq cm 

s distance along surface, cm 
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i' 

t 

u 

V 

O' 

x 

z 

? 

a 

P 

e 

V 
£ 

p 


Pe 



p 

nondimensional surface area, = S/£ 
time, sec 

bubble velocity, cm/sec 
volume of the liquid, cu cm 
nondimensional liquid volume, = V/jo a 
direction parallel to pivot axis, cm 

height above reference horizontal plane, parallel to gravity vec- 
tor, cm 

nondimensional liquid depth; this term- includes cosine of angle be- 
tween axis of rotation and gravity field, z/£ a or a/£ g 

liquid depth in simulator with direction parallel to plates and 
perpendicular to pivot axis (fig. 2), cm 

direction perpendicular to plane of plates (fig. 2), cm 

contact angle, deg or rad 

kinematic viscosity, sq cm/sec 

spacing between plates, cm 

mass density, g/cu cm 

buoyant mass density, - p^,, g/cu cm 
sum over j 

surface tension or surface free energy per unit area, dyne/cm 
angle simulator bed makes from horizontal, rad or deg 

integrals: standard, along curve C, over surface A, over 

volume V 


Subscripts: 
a actual 

B model boundary surface - solid surface not including plates 
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b boundary surface - all solid surfaces enclosing fluid 
f other fluid, in this case air 
i pertaining to interface 

j summing subscript to include all surfaces 

l liquid 

LS liquid- solid interface 

LV liquid-vapor interface 

max maximum 

min minimum 

p surface area of plates within model boundaries 

SLB solid- liquid interface on model boundary 

SV solid-vapor interface 

SVB solid- vapor interface on model boundary 
SVP solid- vapor interface on plates 

s simulated case 
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APPENDIX B 


DISCUSSION OF PLATE SPACING 

The assumption concerning the constancy of the cross-sectional area 
(A(x, a), fig. 2) of the liquid in the planes of the plates requires some addi- 
tional justification. A simple experiment with a drop of wetting liquid (ethyl 
alcohol) "between the two glass plates can he used to show that liquid contacts 
the top and bottom plates when the plate spacing | is below a critical value 
l max - Figure 16 describes the liquid (ethyl alcohol) shape as | increases up 

to i m ax ~ 0.28 centimeter, after which liquid contact with the plates is lost. 
Clearly the variation of the cross-sectional area A(x, a) with p is small for 
large £/|, as long as there is contact. This is assured for £ < l max * A 
nonwetting liquid should give qualitatively similar results; however, the in- 
terface curvature would be reversed. 


The maximum plate spacing £ max that allows liquid to contact the plates 
can be estimated for a wetting liquid by equating the potential energy of the 
liquid before contact is lost to the surface energy of the two free surfaces 
that occur just after contact is lost. Thus 



cos cp)pg Q z dV ^ 2 



for £,/| » 1 


(Bl) 


Since 
uid depth, 
plate (cos 


the cross-sectional area A is essentially independent of the liq- 
the relation simplifies, for constant properties and horizontal 
cp = l), to 


pg 


o 



z dz ^ 2a LV 


(B2) 


or 



(B3) 


Experimental values of £ max are compared in table I with values derived from 
equation (B3). 


TABLE I. - APPROXIMATE VALUES OF | max 
FOR SOME COMMON LIQUIDS 


Liquid 

Experimental, 

Calculated, 


^max' 

^max J 


cm 

cm 

Ethyl alcohol 

0.28 

0.30 

Ethylene glycol 

.33 

.44 

Water 

.43 

.54 


19 



APPENDIX C 


LEVELING 
Parallel Plates 

The velocity of an air bubble injected into the liquid can be used to de- 
termine the angle of parallel plates to the horizontal (fig. 17). The air bub- 
ble should be large enough to contact both plates and should be injected into a 
relatively unrestricted region of the liquid. 

Assume that the predominant forces acting on the nearly cylindrical air 
bubble are the buoyant force and an effective drag force due to form drag of 
the bubble and drag by the plates on the liquid moved by the bubble. The buoy- 
ant force on a bubble of radius a is 

2 

F b = Jta (p z - p f )g D |n sin cp (Cl) 


where 


n - 1 

The effective drag force on the bubble is 


f d = 


jj2 

Vz 2a T £ 


(C2) 


where is the drag coefficient, which is taken to be a function of a 

Reynolds number, of the form 

°D = K <^> e 


(03) 


Assume a modified Reynolds number of the form 


R * = 

R M ~ v l \ a/ 


1/2 


(C4) 


Equating the buoyant and effective drag forces and substituting equations (C3) 
and (C4) into that result give 


U 


2 


tt(p z - p f )s 0 a sin 9 



(C5) 


Experimental data for ethyl alcohol and glycerine at various bed angles and 
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bubble sizes were correlated by equation (C5) when the following values of K 
and e were chosen: K ^ 180 and e ^ - 1. Using the preceding values for K 

and e simplifies equation (C5) to 


U » 


*(p i - p f )g 0 (i) 

" 180 P ' 


3/21 


.V 2 


sin cp 


(C6) 


A plot of experimental data for bed-angle bubble size (i.e., or' sin cp) as a 
function of bubble speed for ethyl alcohol and glycerine (corrected by equa- 
tion (G6) for viscosity) is given in figure 4* 

The accuracy of this method in determining the bed angle at low bed angles, 
as estimated by the scatter of the data in figure 4, is about ±30 percent* At 
low Bond numbers, which correspond to low bed angles, inaccuracies of this order 
in the bed angle would not result in noticeable changes in the liquid shape in 
the two-dimensional model tank* Therefore this method of estimating the low bed 
angles should be adequate for this purpose. 


Plates Canted 


Plate canting (i.e., nonparallel plates) caused by machining inaccuracy in 
the spacers and nonuniform clamping can affect the liquid shape for bed angles 
below about 5x10“ 3 radian* The canting angle error can be compensated for by 
adjusting the level. The combined effect of nonlevel parallel plates and canted 
plates (effective bed angle) can be estimated by the velocity of an air bubble 
by using equation (C6). 
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(b) Side view. 


Figure 3. - Two-dimensional reduced gravity simulator at bed angle of 10' 1 radian (angle exaggerated!. 
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( d) 12. 7-Centi meter-diameter ci rcle at various bed angles. 





Bond number: 
Bed angle: 


460 60 5.15 

2. 29X10* rad 2.97X10' 2 rad 2.56X10" 3 rad 

(e) 7.6-Centimeter-diameter circle with bed angle set so that Bond number is same as for 12.7-centimeter circle. 



Bona number: 
Bed angle: 


169 

0.84X10"* rad 


22 

1.1X10" 2 rad 


1.79 C-65-3608 

8.9X10" 4 rad 




(f) 7.6-Centimeter-diameter circle with same bed angle as 12.7-centimeter circle. 
Figure 5. - Concluded. 














Percent full 
by volume 
(where tank 
is considered 
as surface of 
revolution): 



(a) Simulator results. Bond number (based on diam), approximately 0.02. 





Percent full 
by volume: 

10 


50 


90 C-65-3610 





(b) Drop-tower results. Bond number (based on diam), approximately 0.001. 

Figure 6. - Comparison of liquid shape in partially filled sphere as obtained in drop-tower experiments (ref. 11) to shape observed in 
two-dimensional projection of container in simulator at essentially zero gravity. Sketch is interpretation of adjoining photograph 
by authors of reference 11. 
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Percent full 
by volume 
(where tank 
is considered 
as surface of 
revolution): 



10 


Before flow 
up walls 



Edge effect 


\ Simulator Before flow 

X equilibrium up walls 



50 


(a) Simulator results. 


Bond number (based on tank diam), approximately 0.02. 



Percent full 
by volume: 





(b) Drop-tower results. Bond number (based on diam), approximately 0.006. 


Figure 7. - Comparison of liquid shape in partially filled cylinder as observed in drop-tower experiments (ref. 11) to shape observed in 
two-dimensional projection of container in simulator at essentially zero gravity. 
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as surface of 
revolution): 



10 



(al Simulator results. Bond number (based on length), approximately 0.03. 



Percent full 
by volume: 








C-65-3612 


(b) Drop-tower results. Bond number (based on maximum diam), approximately 0.0015. 


Figure 8. - Comparison of liquid shape in partially filled cone as obtained in drop-tower experiments (ref. 11) to shape observed in 
two-dimensional projection of container in simulator at near zero gravity. 
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Bond number: 467 61.2 4.95 

Bed angle: 0.84X10 -1 rad 1.1X10" 2 rad 0.89X10 -3 rad 

Diameter: 12.7 cm 12.7 cm 12.7 cm 




Bond number: 0.6 

Bed angle: ~10 4 rad 

Diameter: 12.7 cm 



0.02 C -65-3613 

- 10~ 4 rad 
2.54 cm 


Figure 9. - Two-dimensional simulator equilibrium liquid configuration in half-full (based on projected area) 
circular (spherical) tanks at various Bond numbers. 










Bond number: 208 26.2 

Bed angle: 0.84X10" 1 rad l.lX10' 2 rad 

Diameter: 8.45 cm 8.45 cm 



Bond number: -0.123 

Bed angle: -5X10" 5 rad 

Diameter: 8.45 cm 


'-^"Edge effect" 



065-3615 


- 0.022 

- 10"^ rad 
2.54 cm 


2.2 

0.89X10" 3 rad 
8.45 cm 


Figure 11. Two-dimensional simulator equilibrium liquid configuration in half-full (based on projected area) 
rectangular (cylindrical) container (L/D ■ 1.551 at various Bond numbers. 




Bond number: 
Bed angle: 
Tank diameter: 


208 

0.84X10" 1 rad 
8.45 cm 


26.2 

1.1X10" 2 rad 
8.45 cm 



Bond number: 208 

Bed angle: 0.84XKT 1 rad 

Tank diameter: 8.45 cm 



26.2 

1.1X10" 2 rad 
8.45 cm 


Figure 12. - Effect of an internal baffle (cylinder) on simulator i 
rectangular {L/D = 1.55) container at various Bond numbers. 
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2.2 C-65-3616 

0.89X10' 3 rad 


- 0.12 

-5X10' 5 rad 


8.45 cm 


8.45 cm 


ium liquid shape in half-full (based on projected area) 














Bond number: 375 

Bed angle: 0.84X10" 1 rad 

Tank diameter: 12.7 cm 



Bond number.- 3.95 

Bed angle: 0.89X10“ 3 rad 

Tank diameter; 12.7 cm 



47 

1.1X10" 2 rad 
12.7 cm 



— 0. 22 C- 65 -361 7 


~5X10“ 5 rad 
12.7 cm 


(a) Channel (open-end cylinder) axis parallel to direction of gravity field. 

Figure 13. - Two-dimensional simulator equilibrium liquid configuration in 25-percent-fuli 
(based on projected area) circular (spherical) tank at various Bond numbers. 
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Bond number: 375 

Bed angle: 0.84X10' 1 rad 

Tank diameter: 12.7 cm 


47 

1.1 X10" 2 rad 
12.7 cm 



Bond number 
Bed angle 
Tank diameter 


3.95 

0.89X10' 3 rad 
12.7 cm 



-0.22 C-65-3618 


-5X10' 5 rad 
12.7 cm 


(b) Channel (open-end cylinder) axis 45° from direction of gravity field. 


Figure 13. - Continued. 






Bond number: 375 

Bed angle: Q.84X10" 1 rad 

Tank diameter: 12.7 cm 


47 

1.1 X10" 2 rad 
12.7 cm 




Bond number: 
Bed angle: 
Tank diameter: 


3.95 

0.89X10' 3 rad 
12.7 cm 



-0.22 C-65-3619 


-5X10' 5 rad 
12.7 cm 


(c) Channel (open-end cylinder) axis is 90° from direction of gravity field. 


Figure 13. - Concluded. 
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Bond numDer: -0.022 

Bed angle: ~io“ 4 r ad 

Tank diameter: 2.54 cm 


- 0.022 

~ 10" 4 rad 
2.54 cm 

(a) Simulator result. 


- 0.022 
- 10 -4 rad 
2.54 cm 



Bond number: -0.001 

Tank diameter: 5.76 cm 



- 0.001 
5.76 cm 



5.76 cm 

Approaching steady-state zero-gravity 
configuration 


(b) Drop-tower result (ref.l). 

Figure 14. - Further comparison of results of simulator and drop tower. 








S* 0.0254cm 


£ = 0. 152 cm 



Figure 16. - Effect of plate spacing £ on liquid shape between plates for ethyl alcohol 
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